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Abstract. We present soft X-ray spectra of 74 BL Lacer- 
tae objects observed with the PSPC detector on board of 
the ROSAT satellite. The sample contains all BL Lac ob- 
jects detected during the pointed observation phase as a 
target or serendipitously. We have investigated the soft 
X-ray and broad band spectral properties and discuss 
the consequences for the X-ray emission processes. For 
the first time a clear dependence of the X-ray spectral 
steepness on the radio to X-ray spectral energy distribu- 
tion is found: a IX and a x are correlated in the X-ray se- 
lected (XBL) subsample and anticorrelated in the radio 
selected (RBL) subsample. The objects with intermediate 
a rx values thus do have the steepest soft X-ray spectra. 
Simulated PSPC spectra based on a set of simple two 
component multifrequency spectra are in good agreement 
with the measurements and suggest a broad range of syn- 
chrotron cutoff energies. 

We have calculated synchrotron self-Compton beam- 
ing factors for a subsample of radio bright objects and 
find a correlation of the beaming factors Sic with a rx and 
a x . The most extreme RBL objects are very similar to 
flat spectrum radio quasars in all their broad band and 
X-ray properties. 

Key words: Galaxies: active - BL Lacertae objects: gen- 
eral - Radio continuuum: galaxies - X-rays: general 



1. Introduction 

BL Lac objects are active galactic nuclei which by defini- 
tion do not have strong emission lines, are highly variable, 
and show strong polarization in the radio to optical emis- 
sion. Commonly the properties of BL Lacs are explained 
by the concept that the emission in all spectral bands is 
dominated by relativistic jets. Relativistic electrons emit 



synchrotron radiation and may scatter up either syn- 
chrotron photons (synchrotron self-Compton emission) or 
photons from other regions (e.g. from an accretion disc) 
to higher energies. 

Originally most BL Lacs have been found as coun- 
terparts of flat spectrum radio sources (radio selected 
BL Lacs, RBLs), but an increasing number of less radio 
bright objects were discovered by the identification of X- 
ray sources (XBLs). Due to the lack of spectral features 
in the optical no complete optically selected samples ex- 
ist. The total number of BL La c objects in the catalogue 
of Veron-Cetty & Veron ( 1993 ) is less than 200. A new 
catalogue consisting of 233 sources is being published by 
Padovani & Giommi (1996a). The ongoing search for new 
BL Lacs from ROSAT sources (Kock et al. 1996, Nass et 
al. 1996) ) will significantly increase the number of known 
objects. 

Recently BL Lacs have gained great interest due to 
the detection of several objects in high energy 7-rays by 
EGRET on the Compton Gamma Ray Observatory (von 
Montigny et al. 



1995) 



There is some evidence that RBLs and XBLs form dis- 
tinct subclasses as they show a bimodal distribution in 
the plane of the broad band spectral indices a IO versus 
a ox (e.g. Giommi et al. 1990). Regardless of the discovery 
waveband we will use these terms for a distinction of the 
subclasses based on the spectral energy distribution: RBL 
is used for radio bright objects having a rx >0.75 and XBL 
for X-ray bright objects with a rx <0.75. 

There have been several attempts to explain the phys- 
ical differences of the XBL and RBL objects and their re- 
lation to flat spectrum radio quasars ( FSRQs). The con- 
cept that parts of the continuum emission from radio loud 
AGN arises from jets of radiogalaxies more or less aligned 
with the line of sight is widely accepted. 

Based on the study of number count relations and lu- 
minosity functions several authors (e.g. Padovani & Urry 



1990 ) proposed that BL Lacs are the beamed subpopula- 
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tion of FR I galaxies with RBLs having highe r beaming 
factors than XBLs. Ghisellini & Maraschi ( 1989 ) discussed 
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an "accelerating jet" model with lower bulk Lorentz fac- 
tors r in the X-ray emitting regions resulting in broader 
beaming cones for the X-ray emission and narrow radio 



cones. Celotti et al. (1993) developed a "wide jet" model 
with geometrically wider opening angles in the inner, X- 
ray emitting, parts of the jet. Assuming that RBLs have 
smaller viewing angles than XBLs, both models are able 
to explain both the relative numbers and different spectral 
energy distributions of XBLs and RBLs with an intrinsi- 
cally uniform population of objects. Maraschi & Rovetti 
( |l994|) have extended these considerations on FSRQs and 
propose that all radio loud AGN only essentially differ in 
viewing angle and intrinsic power of the central engine. 

An alternative approach to explain the differences be- 
tween XBLs and RBLs was made by Padovani & Giommi 
( |1995|) with a "different energy cutoff" hypothesis. They 
argue that both types form a uniform class of objects 
spanning a wide range in the intrinsic energy distribution 
caused by different cutoff frequencies of the synchrotron 
component. 

The X-ray spectra of RBLs in average were found to 
be significantly steeper than the spectra of (higher red- 
shifted) FSRQs in Einstein IPC (Worrall fc Wilkes |i990|) 



and ROSAT PSPC (Brunner et al. [1994D investigations. 
Furthermore, the X-ray spectral indices of BL Lacs showed 
a broad distribution in both investigations. The mean 
X-ray spectra of XBLs and RBLs were not found to be 
signifi cantly different in Einstein IPC (Worrall & Wilkes 



1990| ), EXOSAT ME + LE (Sam bruna et al. |199J), and 



RO SAT PSPC (Lamer et al. |199J) observations. Ciliegi et 
al. ( 1995 ) found significant steepening of the mean spec- 
trum of XBLs between the soft (0.2-4 keV) and medium 
(2-10 keV) energy X-ray band. 

The ROSAT data archive presently comprises the 
largest X-ray database for BL Lac objects collected by 
a single instrument. In this paper we present the analysis 
of X-ray and broad band spectra of 74 BL Lac objects 
with the X-ray data obtained from pointed ROSAT ob- 
servations. We find a strong interdependence of the X-ray 
spectral index a x with the radio to X-ray energy index 
a rx which we interprete as the signature of two spectral 
components intersecting each other at different frequen- 
cies. 

2. The sample 

Due to the difficulties in the identification and classifica- 
tion of BL Lacs only relatively few objects form complete 
flux limited samples. Nearly complete samples were se- 
lected from radio sources (1 Jy sample, Stickel et al. 1991| ; 
34 objects) and with limited sky cover age fro m soft X-ray 
sources (EMSS sample, Morris et al. 1991 ; 22 objects). 
Our sample comprises all BL Lacs listed in the catalogue of 
Veron-Cetty & Veron ( |l993[ ) and of which ROSAT PSPC 
observations exist in the archives. Only sources which had 
been detected with more than 50 net counts have been 



analysed: 74 objects in total. In case of multiple ROSAT 
observations of an object the longest observation available 
at the time of analysis has been selected. Objects in the 
catalogue which have meanwhile been classified as quasars 
or as radio galaxies have not been included (e.g. Stickel et 



al. 1994). The resulting coverage of various complete sam- 
ples of BL Lac objects is listed in Table |l|. According to 
their radio to X-ray energy distribution 40 objects have 
been assigned to the XBL subsample, 34 objects to the 
RBL subsample (see Sect. 3). 



Table 1. Coverage of complete samples 



sample 


Reference 






observed (total) 


1 Jy (5 GHz) 
EMSS 

S5 (S > 70°) 


Stickel et al. ( 
Morris et al. ( 
Eckart et al. ( 


1991 
1991 
1986 


) 29 (34) 
) 19 (22) 
) 5 (5) 



Throughout this paper the catalogue designations ac- 
cording to Veron-Cetty & Veron (1993) are used. 



3. ROSAT observations and data analysis 

Archival data were taken from the ROSAT data archives 
at MPE (Garching) and at GSFC (Greenbelt). Both the 
author's proprietary data and archival data were reduced 
in the same way using the EXSAS software (Zimmcrmann 
et al. |1993[ ). Tabic | list the objects, ROSAT observation 
request (ROR) numbers, and dates of observations which 
have been analysed. 

The source photons were extracted within a circle of 
radius 100" -200" (depending on the signal to noise ratio) 
and the background determined in an annulus of radii 250" 
and 500" . We produced spectra in the energy range 0.1-2.4 
keV of all objects by binning according to pulse height am- 
plitude, yielding a SNR per spectral bin ranging from 4 for 
the weakest sources and 50 for the strongest. All spectra 
were background subtracted and corrected for telescope 
vignetting, dead time losses, and incomplete extraction of 
source photons with respect to the point spread function. 
The pulse height spectra were then fitted by power law 
spectra combined with the absorption model of Morrison 
and McCammon ( |l983| ). For sources with more than 250 
detected counts both fits with an absorbing column den- 
sity Nn fixed to the galactic value (Stark et al. 1992, Elvis 
et al. 1989| ) and free A% were performed. For the weaker 
sources only fits with fixed 2Vh were obtained. 

In general the latest version of the PSPC detector 
response matrix (nr. 36) has been used for the fits. Except 
for observations carried out before fall 1991, when the gain 
setting of the PSPC was different, an earlier version (nr. 
6) was used. 
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Object 


ROR 


completed 


expos. 


net 


Object 


ROR 


completed 


expos. 


net 




number 


date 


[sec] 


counts 




number 


date 


[sec] 


counts 


PKS 0048-09 


701433 


04-jul-93 


8363 


3368 ± 59 


1207+39W4 


700277 


31-may-91 


38127 


7180 ± 87 


PKS 0118-27 


701427 


10-jul-93 


2635 


366 ±21 


B2 1215+30 


700221 


15-jun-91 


21999 


10141 ± 107 


MS 0122+09 


700976 


16-jul-92 


11173 


208 ± 21 


GB2 1217+34 


700222 


22-may-93 


4281 


282 ± 23 


MS 0158+00 


700972 


24-jul-92 


1646 


1113 ± 35 


1H 1219+301 


701056 


18-jun-92 


10969 


52488 ± 230 


MS 0205+35 


700977 


12-aug-92 


13964 


2644 ± 54 


ON 231 


700223 


14-jun-91 


15223 


3318 ±64 


3C 66A 


600019 


18-aug-91 


24880 


5182 ± 125 


MS 1221+24 


700308 


19-dec-91 


5624 


845 ± 33 


AO 0235+16 


701379 








MS 1229+64 


700974 


03-dec-92 


3034 


2499 ±51 




-...385 


06-aug-93 


14591 


2694 ± 53 


B2 1308+32 


700216 


23-jun-91 


8386 


784 ± 34 


MS 0257+34 


700309 


27-jan-92 


7434 


229 ± 21 


MS 1332-29 


600188 


28-jan-92 


2007 


207 ± 15 


PKS 0301-24 


701503 


10-aug-93 


3412 


936 ±32 


MS 1402+04 


700310 


17-jan-92 


4763 


2876 ± 59 


MS 0317+18 


700971 


24-aug-92 


2917 


412 ± 22 


MS 1407+59 


200406 


10-may-91 


4645 


491 ± 28 


H 0323+02 


700099 


24-aug-91 


25727 


13156 ± 129 


IE 1415+25 


150071 


16-jul-90 


18861 


22242 ± 157 


MS 0331-36 


700921 


05-feb-93 


6214 


424 ± 24 


OQ530 


150046 


19-jul-90 


11516 


2485 ± 55 


1ES 0347-12 


900492 


28-jul-93 


9620 


7817 ±92 


PKS 1424+24 


701504 


ll-jul-93 


3065 


1791 ± 44 


MS 0350-37 


701508 


13-aug-93 


5589 


1227 ± 39 


H 1426+42 


700535 


24-jan-92 


9593 


37670 ± 197 


PKS 0406+12 


700825 


24-feb-93 


15749 


214 ± 16 


MS 1443+63 


700975 


18-jul-92 


4208 


663 ± 29 


1H 0413+00 


701055 


31-aug-92 


9043 


10315 ± 103 


MS 1458+22 


700312 


19-jan-92 


4557 


1694 ±44 


MS 0419+19 


700044 


26-feb-91 


4704 


656 ± 28 


PKS 1519-27 


701430 


17-aug-93 


2552 


262 ± 17 


PKS 0426-38 


701428 


02-aug-93 


4048 


211 ± 18 


4C 14.60 


700834 


28-jan-93 


7215 


282 ± 21 


S5 0454+84 


701058 


17-aug-93 


14850 


63 ± 10 


MS 1552+20 


700973 


26-aug-92 


3562 


2971 ±56 


PKS 0537-44 


700199 


10-apr-91 


2598 


918 ± 31 


MARK 501 


700130 


25-feb-91 


7649 


51 150 ±227 


PKS 0548-32 


700487 


06-mar-92 


9600 


35371 ± 189 


MS 1704+60 


700439 


28-oct-91 


13068 


1090 ± 49 


MS 0607+71 


700111 


13-mar-91 


14086 


247 ± 27 


IZW 187 


701151 


02-oct-92 


4548 


7015 ±85 


S5 0716+71 


700210 


08-mar-91 


21043 


16236 ± 128 


S4 1749+70 


700835 


09-nov-92 


3702 


316±21 


PKS 0735+17 


700829 


28-oct-92 


6684 


592 ±30 


OT081 


700043 


28-mar-92 


9507 


747 ± 31 


MS 0737+74 


800230 


12-mar-93 


8782 


4455 ± 71 


S5 1803+78 


700497 


07-apr-92 


6911 


539 ± 24 


PKS 0820+22 


700831 


29-oct-92 


3728 


64 ±9 


3C 371.0 


700489 


09-apr-92 


10461 


1182± 35 


OJ448 


700832 


04-oct-92 


4821 


59 ± 11 


4C 56.27 


700837 


19-jun-92 


5900 


818 ± 31 


OJ287 


700219 


16-apr-91 


3622 


998 ± 33 


PKS 2005-48 


700488 


27-apr-92 


11487 


33094 ± 177 


B2 0912+29 


700224 


24-apr-91 


2897 


1 172 ± 35 


PKS 2005-48 


701057 


28-oct-92 


11462 


18747 ± 138 


MS 0922+74 


700311 


25-mar-92 


4137 


463 ± 24 


S5 2007+77 


700498 


ll-dec-91 


9525 


443 ± 22 


S4 0954+65 


700042 


17-apr-91 


6772 


343 ± 21 


PKS 2131-02 


701431 


02-nov-93 


7748 


130 ± 15 


GB 1011+49 


700220 


ll-may-91 


5234 


10204 ± 104 


MH 2136-42 


701180 


26-apr-93 


3798 


576 ± 25 


1ES 1028+51 


701544 


22-oct-93 


10783 


23508 ± 157 


MS 2143+07 


700045 


09-may-91 


3526 


748 ± 31 


MARK 421 


700513 


05-may-92 


35838 


5.698 ■ 10" 


PKS 2155-30 


700507 


12-nov-91 


52643 


2.396- 10" 


MARK 180 


700427 


18-nov-91 


3096 


61 17 ± 80 


BL LAC 


700838 


22-dec-92 


2167 


356 ± 19 


PKS 1 144-37 


701434 


07-jul-93 


7747 


852 ± 32 


PKS 2240-26 


701432 


09-nov-93 


3656 


178 ± 16 


B2 1 147+24 


701181 


01-dec-92 


4370 


2842 ± 55 
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We combined the ROSAT measurements with non- 
contemperaneous flux measurements at 5 GHz and in the 
optical V band (both taken from Veron-Cetty & Veron ) 
in order to calculate the broad band spectral indices a rx , 
a ro , and a ox ■ Sources with a rx > 0.75 were assigned to 
the RBL subsample and sources with a xx < 0.75 to the 
XBL subsample. 

Power law spectral indices a are given as energy indices 
(/„ oc v~ a ) throughout the paper. 

We used a maximum likelihood (ML) method to de- 
convolve the measurement errors and the intrinsic distri- 
bution of the X-ray spectral indices and other measured 
parameters when calcu lating mean values and their er- 
rors (see Worrall 1989 for a description of the method 
and Brunner et al. 1994 as a recent application). Assum- 
ing that both the intrinsic distribution of a parameter p 
and the distribution of measurement errors are Gaussian, 
confidence contours of the mean (p) and width ctq of the 
intrinsic distribution of p can be calculated. 

4. Results of spectral analysis 



10 



0.4 



0.6 



1 .0 



R X 



Fig. 1. Distribution of a rx in the full sample (solid line), in the 
EMSS sample (dashed), and in the Uy sample (dotted line). 

The distribution of a rx (Fig. |l|) is double peaked with a 
gap at a rx = 0.6-0.8. As can been seen from Fig. |l], all 
EMSS objects exhibit an XBL energy distribution. Two 
objects of the 1 Jy sample have a rx < 0.75, they therefore 
belong to the XBL sample. This confirms, based on ho- 
mogeneus X-ray data, the relations given by Padovani & 



Giommi ( p95| ). 

Figure |7j shows the interdependence of the calculated 
broad band spectral indices a ox and a ro . Note that the gap 
near a rx = 0.75 is dominated by objects which do neither 
belong to the EMSS sample nor to the 1 Jy sample. This 
indicates that the gap is due to selection effects and may 
be filled in the future when radio or X-ray selected objects 
at lower flux limits will be identified. 



Single power law spectra with photoelectric absorption 
due to the interstellar medium in general yielded accept- 
able fits if the absorbing column density Nn was left free to 
vary. The results of the spectral fits with free Nh and the 
broad band spectral indices are given in Table 3. If a source 
was detected with less than 250 counts, the entry Nn,at is 
omitted and the results with fixed Ah are given. In order 
to investigate whether deviations of the resulting AT H val- 
ues from galactic HI radio measurements are significant we 
calculated the difference AAh for each object. The error of 
the difference was calculated by quadratic addition of the 
X-ray and radio measurement errors using 10 20 cm -2 for 



the Stark et al. ( 1992| ) values and 10 19 cm 2 for the Elvis et 
al. (1989) values. A maximum likelihood analysis of the re- 



sults yields a mean Afe-excess of (0.48 ±0.23) • 10 20 cm~ 2 in 
the XBL sample. In the RBL sample the individual errors 
of the fitted Nn values are generally large and therefore 
no statement about deviations from galactic Ah can be 
made. Note that due to poor energy resolution in the soft 
band of the PSPC an apparent excess of Nn may also be 
caused by a steepening of the intrinsic spectrum. 

As for the fainter objects the statistical error of the 
measured a x is of the same order as the expected varia- 
tions within the sample, we used the maximum likelihood 
method to deconvolve the measurement errors and the in- 
trinsic distribution of the X-ray spectral indices. 

We find (a) = 1.30±0.13, cr G = 0.32±0.I2 for the RBL 
sample (34 objects) and (a) = 1.40±0.09, a G = 0.32±0.07 
for the XBL sample (40 objects). 

Although no significant difference can be found be- 
tween the mean spectra of RBLs and XBLs, there is 
a dependence of the spectral index on the character- 
ising parameter a rx . As can be read from Fig. |^, ob- 
jects with extreme values of a rx on both sides tend to 
have flatter X-ray spectra than intermediate objects. For 
the RBLs (a rx > 0.75) a Spearman Rank test yields an 
anticorrelation of a x and a rx with 99.9% probability. For 
the XBLs (a rx < 0.75) the significance of the positive cor- 
relation is 98%. The steepest X-ray spectra are therefore 
found for objects which fall into the gap between XBLs 
and RBLs (0.6 <a rx < 0.8) for which X-ray spectral in- 
dices up to a x ~2.0 are measured. 

In a previous paper (Brunner et al. 



1994) we stated 



that the mean X-ray spectral index (a x ) of RBLs is similar 
to the mean optical to X-ray spectral slope (a x ), whereas 
the mean X-ray spectrum of FSRQs is significantly flat- 
ter than their optical to X-ray broad band spectrum. 
However, a relatively large dispersion of the differences 
c*ox — o; x was found within the RBL sample. In the larger 
sample presented here a dependence of a ox — a x on a rx is 
visible (Fig. ||) . XBLs show a steepening of the spectrum 
(ttox —Oix < 0), whereas RBLs can exhibit both steepening 
or flattening, depending on a rx . This explains the mean 
(ct ox ~ a x ) being zero with a large dispersion of the in- 
dividual values in an RBL sample. The extreme RBLs 



Table 4. X-ray and broad band spectra 



Object 


JVH,gal 


N H , 6t ± la 


-FlkeV 


tt.l-2.4keV 




6ic 




[10 


!0 cm- 2 ] 


b<Jy] 


±1<T 






PKS 0048 00 


4 92 


a i Q-i_n 53 

1 JTJzW. JJ 


91 


1 70-r-O 1ft 

i . / wzlzw. i w 


82 




PKS 01 1 8-97 


1.60 


9 804-1 50 

Z..OWZT: I .~JJ 


0.19 


i 84-1-0 45 

1 .O^TZtW.^+J 


0.86 




IVliJ U1Z,Z,T^J7 


4 33 




06 


474-0 40 

W.+ 1 ZLZW.+W 


55 




MS 01 58+00 


2 54 


9 ftO-r-0 74 


1 32 


1 33±0 23 


51 




MS 0905+35 

IVliJ \}^\J JTJ J 


6.27 


1 Q QS-r-0 07 

i j . 7 j zlzw. y i 


0.28 


1 504-0 39 

1 . JWZTZW. JZ 


0.54 




3C 66A 


8 89 


1 9 974-0 84 
i z. . z, / nz u.ot 


93 


1 4ft-|-0 39 

1 .+UZIZW. JZ. 


79 


07 


AO 09^5+1 ft 


8 64 


9Q OS-I- 1 77 
Z, J.uoziz I.// 


65 


1 4ftJ-0 40 

1 .+UZtW J-ry 


86 


5 13 


MS 0957+^4 

IVliJ UZ.J ITJt 


10 30 




13 


1 70-1-0 30 
i . / wzlzw. j:? 


64 




PKS 0^01 94 


2 27 


1 434-0 ft1 
i .+ jzlzw.u i 


27 


1 ft8-l-0 31 
i .uoztw. J l 


80 




MS 0^1 7-1-1 8 

IVliJ UJJ / TlO 


10 16 


8 804-9 45 

(I.OU Zt Z . + J 


57 


1 31±0 67 


58 




14 0^9^+09 


8 68" 


Q 30+0 ftO 
o . jwzlzw.uw 


3 21 


1 97-i-O 00 
i .z, / zlzw.w;/ 


57 




MS 0331 3ft 

LVLO UJJJ JU 


1 39 


9 1 S4-1 30 

z, , I j zt I . J w 


34 


1 90-1-0 40 

1 .Z,WZtW. t +:7 


58 




1 FS 0347 1 9 

1 LZ.t5 W J'-t / L Z, 


3 69 


^ 584-0 33 
J. Joztw. J J 


2 56 


1 194-0 00 

1 . 1 Z,ZLZW.W:7 


45 




MS 0350-37 

ivitj uj jw j / 


1.45 


1 59-LO ft1 

1 . JZ,ZtW.U 1 


0.35 


1 17-LO 9ft 
1.1/ _j_yj . ii^yj 


0.61 




PKS 0406+1 9 


14.83 




0.05 


70-1-0 45 

W. / /TU.*tJ 


0.97 




1 H 041 3+00 


9.15" 


Q 1 84-0 SO 
y. 1 oztW. JU 


4 64 


1 554-0 10 

i . j jzitjw. i w 


55 




MS 0410+10 


17.50 


90 i 3-1-99 3 

Z.W. 1 jzczz.j 


0.77 


58-1-0 87 

W . JOZLZW.O / 


0.52 




PKS 049ft 38 
ri\j wm-zu-jo 


1 92 




09 


9 984-0 1 8 
z,.zo znw. i o 


93 




S5 0454+84 


5 60 




02 


874-0 50 
u.o / ztw. jy 


1 01 


3 40 


PKS 0537 44 


4 37 


9 Q44-0 OQ 


81 


i i 84-0 95 


87 


8 97 


PKS 0548 39 
ri\j wjm-o-jz. 


2 49 a 


9 7ft-|-0 14 

Z,. / UZEW. 1 *r 


9 59 


Q54-0 05 


57 




MS 0ft07-i-71 

IVIO UUW / T / 1 


9 63 




08 


1 454-0 48 


70 




S5 0716+71 


3.73 


3.00±0.20 


1.26 


1.75±0.06 


0.77 


1.20 


PKS 0735+17 


4.35" 


2.70±1.32 


0.24 


0.77±0.39 


0.90 


7.70 


MS 0737+74 


3.29 


4.16±0.48 


1.61 


1.39±0.11 


0.54 




PKS 0820+22 


3.92 




0.05 


1.05±0.47 


0.98 




OJ448 


3.89 




0.04 


0.68±0.63 


0.97 




OJ287 


2.75" 


2.52±0.80 


0.62 


1.09±0.29 


0.86 


10.56 


B2 0912+29 


2.82 


1.57±0.55 


0.58 


1.31 ±0.24 


0.72 




MS 0922+74 


1.91 


2.81±1.34 


0.27 


1.03±0.45 


0.53 




IE 0950+49 


1.06 


1.11±0.53 


0.30 


1.45±0.27 


0.52 




S4 0954+65 


4.30 


12.30±3.42 


0.23 


1.06±1.09 


0.89 


5.86 


GB 1011+49 


0.80 


1.35±0.17 


2.15 


1.49±0.08 


0.67 




IE 1028+51 


1.09 


1.33±0.13 


2.55 


1.44±0.05 


0.56 




MARK 421 


1.61" 


2.03±0.01 


272.8 


1.37±0.01 


0.44 


0.29 


MARK 180 


1.27" 


1.87±0.26 


2.62 


1.51±0.10 


0.66 




PKS 1 144-37 


8.74 


10.33±2.87 


0.46 


1.74±0.46 


0.87 





Table 4. - Continued 



Object JVn,gal JVH,fit ± lc -FlkeV a.l_2.4keV «rx SlC 

[10 20 cm- 2 ] [pjy] ±1<7 



B2 1147+245 


1.90 


2.52±0.46 


1.33 


1.25±0.15 


0.76 


0.20 


1207+39W4 


1.82 


3.51±0.34 


0.51 


1.15±0.09 


0.53 




B2 1215+303 


1.60" 


2.16±0.23 


0.49 


1.88±0.08 


0.77 


0.30 


GB2 1217+34 


1.14 


0.94±1.54 


0.10 


0.98±0.71 


0.81 




1H 1219+301 


1.78" 


2.51±0.10 


10.05 


1.22±0.03 


0.48 




ON 231 


1.88 


2.13±0.41 


0.40 


1.24±0.16 


0.81 


0.24 


MS 1221+24 


4.45 


1.76±0.80 


0.20 


1.51 ±0.34 


0.67 




MS 1229+64 


1.96 


3.16±0.56 


2.21 


1.02±0.16 


0.56 




B2 1308+32 


1.10" 


1.31 ±0.77 


0.15 


1.02±0.34 


0.91 


7.72 


MS 1332-29 


3.98 




0.34 


1.14±0.21 


0.59 




MS 1402+04 


2.07 


2.84±0.53 


0.71 


2.03±0.15 


0.57 




MS 1407+59 


1.36 


2.50±1.56 


0.11 


2.07±0.49 


0.67 




IE 1415+25 


2.29 


2.15±0.15 


3.34 


1.19±0.06 


0.55 




OQ530 


1.52 


1.26±0.37 


0.32 


1.09±0.17 


0.85 




PKS 1424+24 


1.65 


5.34±0.81 


1.26 


2.18±0.16 


0.70 




H 1426+42 


2.03 


1.46±0.09 


6.86 


1.05±0.04 


0.49 




MS 1443+63 


1.57 


3.02±1.11 


0.39 


1.09±0.36 


0.58 




MS 1458+22 


3.45 


5.80±0.09 


0.81 


2.26±0.16 


0.59 




PKS 1519-27 


8.72 


10.50±3.06 


0.44 


1.28±0.99 


0.87 


18.69 


4C 14.60 


3.35 


4.02±2.17 


0.12 


0.90±0.72 


0.94 


1.37 




4 06" 




2 62 


ft 74-1-0 1 S 


54 




MARK 501 


1.73 


2.87±0.11 


24.73 


1.77±0.03 


0.62 


1.29 


IE 1704+60 


2.40 


2.50±1.10 


0.16 


1.36±0.35 


0.53 




IZW 187 


2.58" 


3.43±0.34 


3.63 


1.39±0.08 


0.61 


0.01 


S4 1749+70 


3.72 


2.01±1.54 


0.16 


1.23±0.60 


0.89 




B OT081 


9.20 


6.94±4.40 


0.30 


0.47±0.30 


0.88 


14.78 


S5 1803+78 


4.30 


5.16±1.35 


0.23 


1.16±0.34 


0.92 


7.92 


3C 371 


4.67 


3.43±0.81 


0.32 


0.87±0.25 


0.88 


0.81 


4C 56.27 


3.85 


6.46±0.71 


0.52 


0.60±0.34 


0.85 




S5 2007+77 


8.80 


9.50±4.00 


0.20 


0.86±0.77 


0.89 


4.24 


PKS 2005-48 


5.44 


3.84±0.14 


4.00 


1.91±0.03 


0.71 




PKS 2131-02 


3.92 




0.05 


0.85±0.37 


0.99 




MH 2136-42 


2.94 


2.23±1.15 


0.12 


2.20±0.41 


0.77 




MS 2143+07 


4.98 


7.52±1.82 


0.77 


1.52±0.39 


0.63 




PKS 2155-30 


1.35" 


1.41 ±0.04 


46.00 


1.63±0.02 


0.50 




BL LAC 


20.15 


33.84±3.37 


1.46 


2.01 ±1.60 


0.82 


4.00 


PKS 2240-26 


1.54 




0.07 


1.24±0.20 


0.93 





" Elvis etal. 1989, other values: Stark et al. 1992 
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Fig. 2. X-ray spectral index a x versus broad band index a rx , Fig. 3. Model spectra used for the simulations. Dotted: Cs , 
sources with errors Aa x < 0.5 only, filled circles: EGRET de- dashed: Ch , solid: total, shaded area: ROSAT PSPC energy 
tected objects range 



(a rx ~ 0.9) exhibit spectral flattening of the same amount 
as FSRQs (a ox — a x = 0.6, Brunner et al 
The spectral results are available 

3 via WWW (|http://astro.uni 



Table 4. Parameters of model spectra 



199J). 

s a more de- 



tailed version of Tab. 



ix . max 



tuebin gen.de/prepre/), where also an electronic version of 
this paper can be found. 

5. Simulated spectra 

The dependence of the X-ray spectral index on the radio 
to X-ray broad band index ev rx can be explained quali- 
tatively by a two component spectrum as resulting from 



5 • 10 1U 2.00 0.70 0.90 



Sequence of variable parameter a: 

100 200 300 400 500 800 
1000 2000 3000 4000 6000 



synchrotron self-Compton jet models (e.g. Konigl 1981 
Ghiscllini et al 



1985). The correlated variations in a rx and galactic absorption model (Morrison & McCammon |1983[ ) 



cv x are then caused by a varying high frequency cutoff of 
the synchrotron component. 

In order to test this hypothesis we simulated ROSAT 
PSPC spectra with a simple two component spectrum as a 
sum of a parabolically steepening soft component Cs and 
power law hard component Ch (see Fig. |^) : 



Cs(z^) = 1 for v < v\ 

log 10 C s (^ = ^- for v > v x 

V logi a / 

Cyl{v) = N-v- aH 

Below v\ the soft component Cs mimics the flat radio spec- 
trum which BL Lacs have in common; above v\ Cs is a 
parabola in the log v - log /„ plane. The normalization N 
of the hard component was chosen so that Ch (IkeV) and 
Cs (5 GHz) result in a given a rXimax . The extent of the soft 
component can be varied with the parameter a. Figure [3] 
shows the set of calculated spectra using the parameters 
from Tab. [§ 

From the resulting spectra corresponding ROSAT 
PSPC pulse height spectra were determined by applying a 



with Nn = 3 • 10 cm and folding the spectra with the 
PSPC efficiency and detector response matrix. The re- 
sulting pulse height spectra were fitted with an absorbed 
power law model in the same way as the BL Lac spectra. 
Broad band spectral indices a rx , a ro , and cv ox were deter- 
mined from the flux values at 5 GHz, 5517A, and 1 keV 
of each spectrum. 

The locations of simulated and observed spectra in the 
a rx - a x plane are plotted in Fig. ||. Mean values of cv x with 
la errors have been determined in each cv rx interval us- 
ing maximum likelihood contours. We find that the two 
component model is able to reproduce the measured in- 
terdependence of a rx and a x . As curvature of the incident 
photon spectra is able to cause deviations of the JVh values 
resulting from single power law fits, we also compared the 
resulting A^h in the simulated and measured spectra. The 
deviations in A% as derived from the simulations are small 
(< 6 • 10 19 cm -2 ) and thus are hard to detect in individ- 
ual spectra. Averaging the (A^h — A^H.gai) values in bins 
of a rx using the ML method results in an overall excess 
of the measured A^h values (Fig. ||) over the simulations. 
Possibly this excess is caused by intrinsic absorption in 
individual sources. 
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2.5 



2.0 





Fig. 4. Comparison of simulated (filled squares and solid Fig. 6. Flattening (q ox — a x > 0) or steepening (a ox — q x < 0) 
line) and measured (diamonds with error bars) spectra in the of the X-ray spectra with respect to optical to X-ray broad 
a rx - a x plane. The measured spectra were averaged in bins of band spectra vs. a rx . Circles: measured spectra with errors 

A(a ox — a x ) < 0.4), filled squares and solid line: simulated 
spectra. 




1 .0 



o=0.6 - 



Fig. 5. Comparison of simulated (filled squares and solid line) 
and measured (diamonds with error bars) spectra in the a rx - 
ANn plane. The measured AA^h values were averaged in the 
bins indicated by the horizontal bars 



The simulations are able to reproduce the measured 
dependence of the change in the spectral slope between 
the optical to X-ray broad band spectrum and the X-ray 
spectrum, (a ox — ev x ) on cv rx (Fig. ^J). 

Note that the ev rx values of the simulated spectra can- 
not exceed Q rXjmax , which was set to 0.9. Therefore with 
this parameter set the model does not cover objects with 
higher a rx (Figs. |] - ^]). However, the number of objects 
with a rx significantly exceeding 0.9 is small (Fig. |2|). 

In order to test whether the two component models do 
fit the whole radio to X-ray continuum, the broad band 
spectral indices a ox and a IO of simulated and object spec- 
tra were compared (Fig. 0) . The path of the model spectra 
with varying parameter a in the a ox - a ro plane reasonably 
well reproduces the distribution of the measured broad 
band spectra. 



1 1 1 1 1 1 1 1 1 1 

a = 100 
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Fig. 7. Locations of object and simulated broad band spec- 
tra in the a ox - «ro plane. Triangles: EMSS sample, rect- 
angles: 1 Jy sample, circles: others, filled squares and solid 
line: simulated spectra. The division between XBLs and RBLs 
(a rx =0.75) is marked by the dashed line. 



6. Inverse Compton beaming factors 

As inverse Compton beaming factors Sic are often used to 
estimate the viewing angles of radio sources, we investi- 
gated the dependency of Sic on cv rx and a x in our sample. 

The Doppler beaming factor S in a relativistic jet can 
be estimated by the condition that the inverse Compton 
flux from synchrotron self-Compton models must not ex- 
ceed the observed X-ray flux. A method for the calcula- 
tion of Sic from radio brightness temperature and X-ray 
flux was given by Marscher ( 1987 ) : 



Sic = f(a)F B 



Fx 



ln(t/b/V m ) 

(fi+iot a 5+3a 



l/(4+2ct) 



(1 



(1) 
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F m [ Jy] : Synchrotron flux at z^ m [GHz] 

F x [ Jy] : X-ray flux at [keV] 

8d [mas]: VLBI core size 

i>h [Hz] : Synchrotron high frequency cutoff 

a: Optically thin synchrotron spectral index 

f(a) ~ 0.08a + 0.14 (Ghisellini et al. pJ93|) 



We calculated Doppler factors <5ic fro m the compila- 
tion of VLBI data by Ghisellini et al. ( 1993] ) and the 
ROSAT PSPC fluxes at 1 keV; a = 0.75 and v h = 10 14 Hz 
was assumed. The results are given in Table 3. The distri- 
butions of the objects in the a rx - Sic and cv x - Sic planes 
in Fig. H show a strong correlation of Sic with cv rx ■ As the 
subsample with available VLBI data contains predomi- 
nantly RBLs, which show an anticorrelation of a rx and ev x , 
Sic and a x are also anticorrelated. Looking at Eq. (1) two 
reasons may be responsible for the correlation of cn rx and 
Sic ■ 

1. The values of Sic must be considered as lower limits, 
as direct synchrotron emission may contribute to or even 
dominate the X-ray flux and thus the inverse Compton 
flux Fic is overestimated. In this case Sic will be underes- 
timated by the factor 



true /Sic = (F^/Fic) 



-1/5.5 



Assuming that the diversity in 

a rx = (log 10 (F 5G Hz) - log 10 (Fi kcV ))/7.68 

is caused by a more or less energetic synchrotron compo- 
nent, 



Stiuc/S 



ic 



10 



1.40-a r 



results. This function is indicated in Fig. ||a) as solid line. 

2. Both a rx and i5ic depend on the radio flux of the 
source. The total fluxes F 5 qhz determining a rx and the 
VLBI core fluxes F m strongly correlate. Varying the core 
flux _F m and F5GHZ by the same factor while leaving F x 
constant results in 

Sic oc 10 7 68 Q ™ 

as indicated by the dashed line in Fig. || a) . 

Figure || shows that the first possibility, overestima- 
tion of the Compton flux, cannot fully account for the 
measured correlation and viewing angle effects cannot be 
ruled out to cause the differences in Sic ■ 

Ghisellini et al. (1993) noted that BL Lacs show 
a broader distribution of Sic than flat spectrum radio 
quasars ( FSRQs) with a tail towards low values of Sic ■ 
It is apparent from Fig. ^ that one group of objects has 
properties very similar to flat spectrum radio quasars: flat 
X-ray spectra, a rx ~ 0.9, and Sic — 1-10. The correspond- 
ing properties of FS RQs are: (a x ) = 0.59, (a rx ) = 0.88 , 
(Brunner et al. |1994| ), Sic = 1--10 (Ghisellini et al. |1993|) . 
The remaining RBL objects have steeper X-ray spectra, 
lower a rx and lower i5ic • 



10.00 



1 .00 



0.10 



0.01 




X ,ia1 .40arx 

I o 



o ' 



o o 9 



o 
o 



0.4 0.6 0.8 1.0 0.5 1.0 1.5 2.0 

a RX a X 
Fig. 8. a) Inverse Compton beaming factor 5ic vs. a rx • b) 
Sic vs. X-ray spectral index (filled circles: EGRET detected 
objects). 



7. EGRET detected objects 

8 objects of the sample have been detected in high 
energy 7-rays (>100 MeV) by EGRET on the Comp- 
ton Gamma Ray Observatory: AO 0235+16, PKS 0537- 
44, S5 0716+71, S4 0954+65, MARK 421, ON 231, 
P KS 200 5-48 (4-5(7), an d PKS 2155-30 (von Montigny et 
al. |1995| , Vestrand et al. |1995| ). The EGRET objects are 
indicated by filled symbols in Figs. and g. All these 



objects except the low distance XBLs MARK 421 and 
PKS 2155-304 have a rx > 0.7. As can be read from Fig. §, 
the EGRET sources show the same cv rx - a x dependence 
as the remaining objects. Note that not only "FSRQ-likc" 
BL Lac objects have been detected by EGRET, but also 
objects with steep X-ray spectra and moderate Sic ■ 

8. Comparison with other ROSAT studies 

During the refereeing process of this paper we learned 
that a number of ROSAT studies using different BL Lac 
samples are going to be published. In this section we will 
briefly discuss their results in comparison with ours. 

Based on an analysis of 12 objects from the 1 Jy sam- 
ple Comastri et al. ( |1995| ) show that the more radio bright 
objects (a rx > 0.75) in their radio selected sample on aver- 
age have flatter X-ray spectra than the more X-ray bright 
(a rx < 0.75) ones. As their sample covers only the radio 
bright part of the cv rx distribution, this finding fits well to 
our overall picture of the spectrum of BL Lacs, where ob- 
jects of intermediate cv rx do have the steepest X-ray spec- 
tra. 

A detailed investigation of the ROSAT observat ions o f 
the 1 Jy sample has been undertaken by Urry et al. ( 1996 ). 
The above authors both interprete steep X-ray spectra of 
BL Lacs as a sign for synchrotron emission, while flat X- 
ray spectra should be dominated by self-Compton emis- 
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sion. Perlman et al. (1996) performed a ROSAT investi- 
gation of the EMSS XBL sample and found a distribution 
of X-ray spectra similar to the 1 Jy sample. 

By considering the whole sample of BL Lac objects 
we are able to verify the view, that BL Lacs except the 
extreme RBLs are dominated by synchrotron emission. 
We show that the extreme XBLs have flat X-ray spectra 
caused by synchrotron spectra with cutoff energies beyond 
the soft X-ray band. Our spectral simulations yield a good 
measure for the synchrotron cutoff energy and show that 
the range of cutoff energies is large. 

Padovani&Giommi (1996b) have calculated X-ray 



spectral indices in a large sample of BL Lacs from hardness 
ratios prov ided by the ROSAT WGA catalogue (White 
et al. 1994) and with this method obtain a similar depen- 
dency of a x on the radio to X-ray flux ratio as we do. 

9. Discussion 

We find that the broad distribution of spectral slopes in 
the soft X-ray spectra of BL Lacs is due to a strong de- 
pendence of a x on the broad band spectral index a rx . Ob- 
jects with extreme values of cv rx exhibit flat X-ray spec- 
tra, whereas intermediate objects have steeper spectra. 
The symmetry of this dependence prevented the detection 
of significant differences between the mean X-ray spectra 
of XBLs and RBLs in previous investigations (Worrall & 
Sambruna et al. 1 

with simulated 



1994) 



Wilkes 1990, Sambruna et al. 1994, Lamer et al. 

By comparison with simulated PSPC spectra we 
showed that a two component model with a hard power 
law (a — 0.7) component and a steepening soft compo- 
nent is appropriate to explain the observed spectra. The 
frequency where the components intersect each other is 
below the soft X-ray band for extreme RBLs and crosses 
the energy band of the ROSAT PSPC with declining 
cv rx . In the framework of the SSC models this means that 
Compton emission causes the flat X-ray spectra of extreme 
RBLs, whereas the likewise flat X-ray spectra of extreme 
XBLs are due to synchrotron emission. The steep X-ray 
spectra of objects with intermediate spectral energy distri- 
bution (0.6 <a rx < 0.8) represent the first direct evidence 
of the synchrotron high energy cutoff. 



Padovani & Giommi (1995) explain the different spec- 
tral energy distributions of XBLs and RBLs by different 
energy cutoffs of the synchrotron spectra. They postulate 
the cutoff energy being intrinsic properties of the sources 
without discussing the physics of the emission processes. 
This scenario also is the most straightforward explana- 
tion for our findings, including the correlation of a rx and 
ev x for the XBL subsample. The wide range in synchrotron 
cutoff energies, and consequently the cutoff in the energy 
spectrum of the relativistic electrons, has to be explained. 
The SSC cooling of the jet electrons may be more efficient 
in the more powerful jets of FSRQs an d RBL s than in 
the jets of XBLs. Ghisellini & Maraschi ( |1994|) proposed 
a more rapid cooling of jet electrons in FSRQs by exter- 



nal UV photons compared to jets of BL Lac objects. It is 
conceivable that RBLs are intermediate objects between 
XBLs and FSRQs regarding the ambient photon density. 

The more physically motivated beaming models, such 
as the "accelerating jet" model (Ghisellini & Maraschi 



1989) and the "wide jet" model (Celotti et al. 1993) do 



not as naturally satisfy our data. The crossover frequen- 
cies of soft and hard components in the spectra calculated 



by Ghisellini and Maraschi (1989) do not span a sufficient 



range and do not move across the soft X-ray range when 
tilting the viewing angle, as required by the new data. 
Nevertheless, further tuning of the free parameters may 
provide spectra which are in accordance with the mea- 
surements. The ROSAT spectra therefore are suitable to 
constrain the beaming models. 
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